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ABSTRACT 



We study the effect of high magnetic fields on the particle population of 
neutron stars using an extended hadronic and quark SU(3) non-linear realization 
of the sigma model. In this model the degrees of freedom change naturally from 
hadrons to quarks as the density and/or temperature increases. The effects of 
high magnetic fields and anomalous magnetic moment (AMM) are visible in the 
macroscopic properties of the star, such as mass, adiabatic index, moment of 
inertia, and cooling curves. 

Subject headings: stars: neutron — stars: magnetars — stars: interiors — stars: 
magnetic field — stars: evolution 



-3 - 



Introduction 



Magnetars are compact stars th at have surface magnetic fields up to 10 14 — 10 15 G 



( jPaczynski 



1992 



Ibrahim et al 



20041 ). Such fields can be estimated from observations of the 
star's period and period derivative. According to Virial theorem estimates, neutron stars 
could have a central magnetic field as large as 10 18 or 10 19 G in their interiors. The precise 
calculation of this limit is complicated, since all energies to be weighed with the magnetic 
energy (the one from matter and the gravitational one) also depend on the magnetic 
field. For this reason the limit is different for different equations of state. Furthermore, 
due to the asymmetry introduced by the magnetic field in the z-direction, the pressure 
becomes different in the parallel and perpendicular directions. In reality, depending on 
the magnitude of the magnetic field, the parallel pressure becomes muc h smaller than the 



perpendicular one, and in extreme cases, can go locally to zero. In ref. ( jFerrer et al. 



2010), 



as well as in our calculations, this limit was found to be around 10 18 - 10 19 G. Beyond this 
value strong instabilities can occur, which indicates that further investigation is necessary. 



References ( IBocquet et al.l Il995l ; ICardall et al.l 120011 ) solve the Einstein and Maxwell 
equations self-consistently for non-rotating and rotating stars. They study axisymmetric 
and poloidal magnetic field configurations and find B ~ 0.1 to 4.2 x 10 18 G as limits for the 
star central magnetic field. The range comes from the use of different hadronic equations of 
state, from a simple politropic to a hyperonic relativistic one. The equations of state used 
do not take into account magnetic field modifications, but the formalism provides reliable 
results for star masses and radii because it takes into consideratio n different pressures in 



2002|) 



different directions of the star. It is important to notice that Ref (IBroderick et al. 
states that the use of different current functions and symmetries different that the axial one 
together with effects of the magnetic field on the equation of state can alter these limits. 



On the other hand, it was shown by Ref. (jKurkela et al. 



2010 



bj) that using well 
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accepted hadronic EOSs together with one based on first principles QCD calculations, there 
is indeed a phase transition to deconfined matter already at a few times saturation density 
inside neutron stars. Such a phase transition can be a sharp first order transition or include 
a mixed phase, depending on the surface tension between the phases. In this work we 
study the effect of the latter. One important question we will try to address is if the effects 
of strong magnetic fields in the deconfinem ent phase transition are strong enough to be 



observed. Other studies along this line are (jChakrabarty 



Rabhi et al 



2009 



Paulucci et al 



1996 



Bandyopadhyay et al. 



1997 



201ll ). In order to do this in a realistic way we adopt as 



our microscopic model a chiral approach that includes hadronic as well as quark de 



grees of 



20101). Here, 



freedom in a unified description that was developed in (iDexheimer fc Schramm! : 
besides the analysis of the mass-radius diagram (in one dimension only), we will extend our 
study to the adiabatic index and moment of inertia of the stars. Furthermore, we will, for 
the first time, investigate the thermal evolution of highly magnetized hybrid stars, which 
will allow us to assess the effect of a high magnetic field on the cooling of the star. 



The Model 



Chiral sigma models are effective quantum relativistic models that describe hadrons 
interacting via meson exchange and, most importantly, are constructed from symmetry 
relations. They are constructed in a chirally invariant manner as the particle masses originate 
from interactions with the medium and, therefore, go to zero at high density/temperature. 



the sigma mode 



gives 



1998 



Adopting the nonlinear realization o 
the widely used linear sigma model (jPapazoglou et al. 
in be tter agreement with nuclear physics results ( jPapazoglou et al. 
20091 ). 



a significant improvem ent to 



Lcnagh an et al. 



1999 



2000) as it is 



Bonanno fc Drago 



The Lagrangian density of the SU(3) non-linear sigma model in the mean field approx- 
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imation constrained further by astrophysics data can be found in (IDexheimer fe Schramm 



20081 ). A recent extension of this mod el also includes quarks as dynamical degrees of 



freedom (IDexheimer &: Schramm 



2010). In this version, the degrees of freedom change 
due to the effective masses of the baryons and quarks. Their masses are generated by the 
scalar mesons (a, isovector 5, strange (), except for a small explicit mass term Mq and the 
term containing t he field $, related to deconfinement in analogy with the Polyakov loop 



(IFukushima 



20041 ). The baryon and quark masses are given by 



m *b = 9BaCr + g B 6T 3 5 + g B cC + M o B + 9b<s>® 2 



M* 



9qaO- + g q sr 3 5 + g q( ( + M 0q + g qq> (\ - $), 



(2) 



2010). 



where the values for the coupling constants can be found in (IDexheimer fc Schramm 
With the increase of density/temperature, the a field (non-strange chiral condensate) 
decreases its value, causing the effective masses of the particles to decrease towards 
chiral symmetry restoration. The field $ assumes non-zero values with the increase of 
temperature/density and, due to its presence in the baryons effective mass (Eq. (JT])), 
suppresses their presence. On the other hand, the presence of the $ field in the effective 
mass of the quarks, included with a negative sign (Eq. ([2])), ensures that they will not be 
present at low temperatures/densities. 

Both transitions, chiral symmetry restoration and deconfinement happen at the same 
density/temperature. Such a fact comes from the correlation of these quantities in the 
effective masses of the particles. The potential for $ reads 



U = (a T 4 + a lA 4 + ojTX)^ + a 3 T 4 log (1 - 6$ 2 + 8<^ - 3$ 4 



(3) 



and allows the model to be used to study low, as well as high temperature environments, 
and make p redictions for high-energy heay y-ion physics and comparisons with lattice QCD 



results (see (IDexheimer fc Schramm 



201(J) for details). 
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3. Inclusion of Magnetic Field 



In order to have a more detailed description of neutron stars, we include magnetic fields 
that, while being in the z-direction, are not constant. The effective magnetic dipole field 
B* increases with baryonic chemical potential going from a surface value of B sur f = 69.25 
MeV 2 = 10 15 G to a central value B c 



B* = B 



surf 



BJl 



o"\ 939 )\ 



(4) 



where /is is understood to be given in MeV. In this way, there is no discontinuity in the 
magnetic field or in the increase of the magnetic field at the phase transition. Such an 
unphysical discontinuity would be present if we had chosen the effective magnetic field to 
be a function of baryonic density, as is usually done. The constant b is chosen to reproduce 



(in the absence of guar 



(IBandyopadhyay et al 



ks) t 



i e effective magn etic field curve as a function of density from 



1997 



Mao et al. 



20031 ). As can be seen in Fig [TJ even with the 



use of extremely high central magnetic fields (recall that using Gaussian natural units 1 
MeV 2 = 1.44 x 10 13 G), the values for the effective magnetic fields only become extreme 
at very high baryonic chemical potentials that correspond to the inner core of the star. 
Such a depe ndence prevents effects like the spontaneous creation of electron-positron pairs 



in the star (jSchwinger 



195ll ). Still, the highest central magnetic field considered in this 



work is considerably above the limit est ablished by hydrostatic stability from references 



( IBocquet et al 



1995 



Cardall et al. 



20011 ) and is only shown to illustrate the influence of an 



extreme high magnetic field in the EOS. 

The magnetic field in the z-direction forces the eigenstates in the x and y directions of 
the charged particles to be quantized into Landau levels v 



E* 



kl 



mf + 2u\ qi \B* - s^B* 



(5) 



where hi is the fermi momentum, qi the charge and Sj the spin of each baryon or quark. 
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The last term comes from the anomalous magnetic moment (AMM) of the particle that 
splits the energy levels with respect to the alignment/anti-alignment of the spin with the 
magnetic field. The AMM also modifies the energy levels of the uncharged particles 



E* 



rn 



* 2 



iKiB*) . 



The constants Ki determine the tensorial coupling strength of baryons with the 



electromagnetic field tensor and have values k„ = 1.79, K n = —1.91, k,\ 



-0.61, k+ 



(6) 



1.67, 



NX 



1.61, k e = —0.38, k| = —1.25, k h = 0.06. The sign of Ki determines the preferred 



orientation of the spin with the magnetic field. The sum over the Landau levels v runs up 
to a maximum value, beyond which the momentum of the particles in the z-direction would 
be imaginary 



E* s 2 + s iKi B* 



rn 



2 



2\qi\B* 

We choose to include in our calculations the AMM effect for the hadrons only, since the 
coupling strength of the particles k$ depends on the corresponding magnetic moment, 
that up to now is not fully understood for the quarks. Furthermore, it is stated in Ref. 



(7) 



([Weinberg 



moment. 


Jbr c 


Suh et al. 


2001 



1990|) , that quarks in the constituent quark model have no anomalous magnet ic 



'or calculations including AMM effects for the qua rks see ( iChakrabarty 



1996 



Perez Martinez et al. 



2005 



Felipe et al. 



20081 ). The AMM for the electrons 



is also not taken into account as its effect is negligibly small. The AMM removes the spin 
degeneracy of the particles and their energy levels are further split in two levels each, which 
increases even more the particle chemical potentials. 

As can be seen in Fig. [2j a phase transition (chiral symmetry restoration and 
deconfinement) which is of first order at zero or small temperatures is delayed with the 
inclusion of high magnetic field. This effect is due to a stronger stiffenening of the hadronic 
part of the EOS (as the magnetic field increases). It is important to notice that at very low 
densities, the hadronic EOS of non-interacting matter gets softer in the presence of high 
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magnetic field, in agreement with ( IBroderick et al.ll2000l . 120021 ) and many other papers. In 
our case this effect is smaller because the magnetic field is only large at high densities. The 
EOS is shown in Fig. [3] for the case in which global charge neutrality is allowed and a mixed 
phase appears. We do not show the pure quark phase part, since in this model it is present 
only at very high densities that cannot be reached in the interior of the neutron star. 

The delay of the phase transition is further increased when the anomalous magnetic 
moment is considered, as it only affects the hadronic phase rendering its EOS stiffer (due 
an increase in the chemical potential of the baryons). If global charge neutrality is assumed, 
the location of the mixed phase does not change substantially, except for the highest 
magnetic field considered with AMM, when the entire mixed phase is pushed to higher 
chemical potentials. 

The same effect of confinement and/or chiral symmetry enhancement i n the presence 



2009 


b; 


Avancini et al. 


201: 


.) this was shown for high 


In references ( 


D'Elia et al. 


2010; 


Mizher et al. 


2010; 



i density and small or zero temperature. 



Avancini et al. 



2011 



Gatto &; Ruggieri 



20 111 ) this was shown for small or zero density and high temperature. Together, these 
features show the importance of the study of the influence of high magnetic fields and 
deconfinement /chiral symmetry restoration in compact stars as a part of a greater picture 
that forms the whole QCD phase diagram and contains heavy-ion experiments in the other 
extreme. 

Fig. H] shows the baryonic density of fermions (quark number densities are divided by 3). 
In the mixed phase the hadrons disappear smoothly and the quarks appear smoothly. The 
hyperons, in spite of being included in the calculation, are suppressed by the appearance 
of the quark phase. Only a very small amount of A baryons appear right before the phase 
transition. The density of electrons and muons is significant in the hadronic phase but not 
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in the quark phase. The reason for this behavior is that, because the down and strange 
quarks are also negatively charged, there is no necessity for the presence of electrons to 
maintain charge neutrality, and only a small amount of leptons remains to ensure beta 
equilibrium. The strange quarks appear after the other quarks and do not make substantial 
changes in the system. 

Fig. [5] shows the change in population when a central magnetic field of 5 x 10 5 MeV 2 
with AMM is considered. The wiggles in the charge particle densities come from the Landau 
levels, more precisely when the Fermi energy of the particles cross the discrete threshold of 
a Landau level. The charged particles are enhanced (as their chemical potentials increase 
with B), which is especially visible in the amount of electrons in the quark phase. The 
remaining hyperons are even more suppressed due to the increase in the proton density. All 
these effects are further enhanced when higher central magnetic fields are considered, like 
B c = lx 10 6 MeV 2 . 



4. Macroscopic Properties 

A key point of this investigation is to determine whether the changes due to the 
presence of magnetic field in the EOS are strong enough to affect observable properties of 
the stars. To answer this question, we are going to analyze the adiabatic index, moment of 
inertia, mass-radius diagram and thermal evolution for different magnetic fields. 



4.1. Adiabatic Index 

We begin with the adiabatic index (Fig. [6]). The discontinuities in the B = curve 
show the appearance of the A's and the beginning of the mixed phase. The extra peaks 
in the finite magnetic field with AMM curves show the Landau level thresholds that also 
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can be seen in the population plot. Note that the Landau level t hresholds appear in the 



201 lh . these 



hadronic phase as well as in the mixed phase. As pointed out in ( IRyu et al. 
rapid changes in the pressure can cause instabilities in the star that might cause star-quakes, 
glitches and giant flares. 



4.2. Moment of Inertia 

In Fig. [7] we show the moment of inertia as function of radius for several of the 
investigated stellar sequences taking into account the AMM. We see that for the maximum 
central magnetic field considered (1 x 10 6 MeV 2 ), there is an increase of about 10% from the 
unmagnetized case (for the maximum mass star). For B c = 5 x 10 5 MeV 2 this increase is of 
about ~ 2%. This result indicates that if the magnetic field of the object is changing, the 
moment of inertia will be modified, which might affect the spin evolution of the object, and 
mi ght allow one to make a connection between AXP emission and glitches (see discussion 



in ( Ryu et al. 



201 if )). For lower central magnetic fields, the moments of inertia are almost 



the same as in the non-magnetized case. 



4.3. Mass radius diagram 



As a first approach to the problem, the possible neutron star masses and radii are 



calculated b y solving the 



static stars (iTolman 



1939 



blman-Oppenheimer-Vo 



Oppenheimer &: Volkofi 



koff equations for spherical isotropical 



19391 ). In Fig. [HI be sides our equation of 



19711). 



state for the core, a separate equation of state was used for the crust (iBaym et al 
The maximum mass supported against gravity is higher for higher magnetic fields and 
even higher when the AMM is included. It is important to note that in this model pure 



quark matter only appears at very high densities, that correspond to the unstable branch 
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of the mass-radius diagram. Thus, only hadronic and "mixed" matter appear in the star. 
Similar results for stars containing only hadronic and "mixed matter" were also found 
in a calculation using the Brueckner-Hartr ee-Fock model fo r the hadronic phase and the 



Dyson-Schwinger approach for the quarks (IChen et al. 



2011 



So far, the magnetic field pressure was not considered and only the influence of the 
magnetic field on the energy levels of the particles was taken into account. The problem in 
including the magnetic pressure is that it has different values in the directions parallel and 
perpendicular to the field, as can be seen in the electromagnetic energy-momentum tensor 

- u u u u ^ 

1 



rpfll/ 

_ 4:71 



V 

















1 r>*2 
2 














lg*2 

















J 



where the first term is related to the energy density and the other three to the pressure 
i n the x, y and z directions. This problem was pointed out in many papers such as 



Chaichian et al 



2000 



Perez Martinez et al 



2003 



'erez Martinez et al 



2005 



Felipe et al 



2008 



Perez Martinez et al. 



2008 



Ferrer et al 



2010 



Orsaria et al. 



2011 



Paulucci et al 



201ll ). In these papers, there is also a term coming from the magnetic dipole interaction, 



which is linear in B and is therefore subleading at high magnetic fields. 

A consistent inclusion of the macroscopic magnetic pressure requires a 2D calculat ion. 



As mentioned before, such calculations exist ( iBocquet et al. 



1995 



Cardall et al 



200lh : 



however, most calculations of the macroscopic properties of magnetized stars assume 
isotropy and consider the pressure term to be positive in all d i rections or negative in a 



direc tions. It was pointed out by (IBandyopadhyay et al. 



1997 



Sinha fc Mukhopadhyay 



2010) that negative pressure in all directions constrains the maximum values that can 
be used for the magnetic field to lower values. In our case, B c = 5 x 10 5 MeV 2 would 
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already cause such in stability at very high densities. A third option, suggested by 



( IBednarek et al 



20031 ) uses, as monopole approximation of the energy-momentum tensor, 
the average between the three pressures and adds +B 2 /2Atc for the magnetic pressure in all 
directions. We show the difference caused by either choosing the pressure to be positive 
in all directions or using the average pressure in the mass-radius diagram. Obviously the 
first case is unphysical as the magnetic pressure dominates over any other contribution. We 
speculate that the correct value would be close to the third option that uses the average 
of the pressures in different directions, but any mor e precise statement r equires the use of 



20111 ). Work along this 



axisymmetric equations, as was already pointed by (IPaulucci et al.l 
line using our hybrid matter EOS is ongoing and will allow for a more exact estimate of the 
increase of the star mass based on the magnitude of the magnetic field of the magnetar. 



4.4. Thermal Evolution 

We now turn our attention to the thermal evolution of hybrid stars, whose microscopic 
composition is given by the model described in this paper. The cooling of compact stars is 



given by the thermal balance and thermal energy transport equation (G = c = 1) (IWeber 



1999|) 



d(le») 1 ^ , „ d(Te 



dm pa/1 - 2m/r \ dt 

d(Te^>) (Ze* 



ce^ + c-Ss-^ ) > (9) 

(10) 



dm 16iT 2 r 4 K,p^l - 2m/r ' 

In Eqs. (191 — (TTOl) the structure of the star is given by the variables r, p(r), and m(r), 
that represent the radial distance, the energy density, the metric function, and the stellar 
mass, respectively. The thermal variables are given by the temperature T(r,t), luminosity 
l(r,t), neutrino emissivity e v (r,T), thermal conductivity n(r,T), and specific heat c v (r,T). 
The boundary conditions of (Q and fflOl) are determined by the luminosity at the stellar 
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center and at the surface. The luminosity vanishes at the stellar center since there is no 
heat flux there. At the surface, the luminosity is defined by t he relationship between th e 



mant l e temperature and the temperature o utside of the star (IGudmundsson et al. 



1983 



Page et al. 



2006 



Blaschke et al.l 



1982 



20001 ) . For the hadronic phase, we have considered 
the following neutrino emission processes: direct Urea, modified Urea and bremsstrahlung 
processes; whereas for the quark phase, the processes taken into account are the quark 



direct Urea, quark modified Urea, and quark bremsstrahl ung processes. Petal 



emiss ivities of such processes can be found in references (jYakovlev et al. 



2001 



s about the 



Iwamoto 



982 ). The specifi c heat of the hadrons is the usual specific heat of fermions, as described in 



(IPage et al. 



in ( jlwamoto 



2004). As for the quarks, we use the expression for the specific heat calculated 



19821 ). The t hermal conductivity 



calculated, respectively in ([Flowers fc Itoh 



1981 



or the hadroni c and quark matter were 



Haensel 



19911 ). and in this work we follow 



the results presented in these references. 

We have calculated the cooling for magnetized hybrid stars (with anomalous magnetic 
moment of the hadrons included) for three different masses: 1.1, 1.4 and 1.93 M . For each 
hybrid star mass we considered three values for the magnetic field (1 x 10 5 , 5 x 10 5 , 1 x 10 6 
MeV 2 ). In addition to that we have also considered the effect of including the magnetic 
field pressure in the EoS. 

The cooling of 1.1 M & stars is shown in Fig. [9j One can see that in this case the star 
exhibits a slow cooling, which agrees relatively well with the observed data. Furthermore, 
we can also conclude that in this case the magnetic field has no substantial effect on the 
thermal evolution of the object, and neither does the inclusion of the magnetic pressure 
term in the EoS. 

The situation is similar for stars with higher masses, as can be seen in Figs. [10] and [HI 
which show the thermal evolution of stars with 1.4 and 1.93 M Q , respectively. We see that 
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in this case the modifications introduced by the magnetic field in the composition are not 
enough to alter the cooling (as was the case for lower mass stars). This picture changes, 
however, if one introduces the effect of the magnetic pressure. As shown in Figs. [TUland[TT| 
the extra _B* 2 /87r leads to slow cooling of a star that would otherwise exhibit fast cooling. 
This stems from the fact that adding an extra source of pressure stiffens the EoS. Hence, 
stars with higher masses possess smaller central densities and, therefore, smaller proton 
fractions. The small proton fraction will hinder the otherwise present direct Urea process, 
thus leading to a slow cooling. The results for the inclusion of the magnetic pressure of 
B* 2 /24n are qualitatively the same. In this case, however, the stiffening of the EoS is more 
mild, and thus the effect is only appreciable for moderate masses, as seen in Fig. [TUl where 
the thermal evolution is the same as for the magnetic pressure of B* 2 /8tt. For the 1.93 M & 
star we see that the stiffening of the EoS is not enough to hinder the direct Urea process 
completely, and cooling of neutron stars with magnetic pressure of B* 2 /24tt is only slightly 
slower. 



5. Conclusions 

To model magnetars, we use an effective model that includes hadronic and quarks 
degrees of freedom. As the density increases, the order parameters signal the deconfinement 
and chiral symmetry phase transitions, which happen more towards the center of the star 
for higher magnetic fields. The stiffness of the EOS, and the consequent star masses, depend 
on the chosen central magnetic field. It is clear that a higher B c allows for more massive 
stars, but the quantitative analysis of how more massive the star can become requires 
the use of a 2D solution of Einstein's equations which takes into account the breaking 
of spherical symmetry by the magnetic field. Such formalism will also allow to calculate 
the exact hydrodynamic limit for the magnetic field for our model. Efforts along this line 
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using our EOS are currently being pursued. In addition to the effect on the mass-radius 
relationship, in this work we have shown that the presence of high magnetic fields can result 
in peaks in the adiabatic index and can modify the star's moment of inertia. 

In our work we introduce an expression for the effective magnetic field B* that 
increases with baryonic chemical potential, and therefore avoids a discontinuity in the phase 
transition region. This allows B* to increase smoothly from a surface value up to a chosen 
central one (B c ). Furthermore, our investigation indicates that the composition changes, 
introduced by the magnetic field, are not high enough to alter the thermal evolution of 
hybrid stars. We have found, however, that by including the magnetic pressure (_B* 2 /87r), 
the resulting stiffening of the EoS allows for stars with higher masses (1.4-1.93 M©) leading 
to a slower cooling, which is in better agreement with the observed data. This is an 
important result, which indicates the importance of the introduction of the macroscopic 
pressure correctly into the structure of the star. We are cautious in interpreting these 
results, since the appropriate treatment of both the structure, and the thermal evolution 
of stars with such high magnetic fields requires a full two-dimensional study. We will save 
this investigation for a future work, as it seems clear to us that the effect of the magnetic 
pressure on the structure and thermal evolution are substantial. Although calculations 
which address this problem using model equations of state have already been performed, we 
believe that a well-motivated microscopic equation of state that also includes the effects of 
the star's magnetic field is important for fully understanding the behavior of magnetars. 

Subsequently, we would like to expand our study to finite temperature, in order to 
have a complete phase diagram (temperature as a function of baryonic chemical potential) 
for high magnetic field. A phase diagram for T = was al ready constructed using the 



extended non-linear realization of the SU(3) sigma model in (iDexheimer fc Schramm 



2010l ). 



Such an analysis is important, since it connects the physics at high temperature and small 
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densitiessuch as those created in heavy ion collisions, with the p hysics inside neutro n 



19961 ) but 



stars. Such a finite-!? phase diagram was already constructed by (iChakrabartyl 
we believe that we can give further insight into the matter with our model, since it also 
includes chiral symmetry restoration and allows for a smooth chiral and deconfinement 
crossover transition at small densities. 



We are grateful to Constanca Providencia, Milva Orsaria, Fridolin Weber, Debora 
Menezes and Aurora Perez Martinez for the fruitful discussions on the subject of strong 
magnetic fields. R. N. acknowledges financial support from the LOEWE program HIC for 
FAIR. 



-17- 
REFERENCES 

Avancini, S. S., Menezes, D. P., & Providencia, C. 2011, Phys. Rev. C, 83, 065805 
Bandyopadhyay, D., Chakrabarty, S., & Pal, S. 1997, Phys. Rev. Lett., 79, 2176 
Baym, G., Pethick, C, & Sutherland, P. 1971, Astrophys.J., 170, 299 

Bednarek, I., Brzezina, A., Manka, R., & Zastawny-Kubica, M. 2003, Nucl.Phys., A716, 245 
Blaschke, D., Klahn, T., & Voskresensky, D. 2000, The Astrophysical Journal, 533, 406412 
Bocquet, M., Bonazzola, S., Gourgoulhon, E., & Novak, J. 1995, A&A, 301, 757 
Bonanno, L., & Drago, A. 2009, Phys. Rev., C79, 045801 
Boomsma, J. K., & Boer, D. 2009, Phys. Rev. D, 80, 034019 
Broderick, A., Prakash, M., & Lattimer, J. M. 2000 

Broderick, A. E., Prakash, M., & Lattimer, J. M. 2002, Phys. Lett., B531, 167 
Cardall, C. Y., Prakash, M., & Lattimer, J. M. 2001, ApJ, 554, 322 

Chaichian, M., Masood, S. S., Montonen, C, Perez Martinez, A., & Perez Rojas, H. 2000, 
Physical Review Letters, 84, 5261 

Chakrabarty, S. 1996, Phys. Rev., D54, 1306 

Chen, H., Baldo, M., Burgio, G. F., & Schulze, H. -. 2011. larXiv: 1107. 24971 
D'Elia, M., Mukherjee, S., & Sanfilippo, F. 2010, Phys. Rev. D, 82, 051501 
Dexheimer, V., & Schramm, S. 2008, Astrophys. J., 683, 943 
Dexheimer, V. A., & Schramm, S. 2010, Phys. Rev., C81, 045201 



-18 - 

Felipe, R. G., Martinez, A. P., Rojas, H. P., & Orsaria, M. 2008, Phys. Rev., C77, 015807 

Ferrer, E. J., de la Incera, V., Keith, J. P., Portillo, L, & Springsteen, P. P. 2010, Phys. 
Rev., C82, 065802 

Flowers, E., & Itoh, N. 1981, The Astrophysical Journal, 250, 750 

Fukushima, K. 2004, Phys. Lett., B591, 277 

Gatto, R., & Ruggieri, M. 2011, Phys. Rev. D, 83, 034016 

Gudmundsson, E. H., Pethick, C. J., & Epstein, R. I. 1982, The Astrophysical Journal, 259, 
L19 

Gudmundsson, E. H., Pethick, C. J., & Epstein, R. I. 1983, The Astrophysical Journal, 272, 
286 

Haensel, P. 1991, Nuclear Physics B - Proceedings Supplements, 24, 23 
Ibrahim, A. I., et al. 2004, AIP Conf. Proa, 714, 294 
Iwamoto, N. 1982, Annals of Physics, 141, 1 

Kurkela, A., Romatschke, P., & Vuorinen, A. 2010, Phys. Rev. D, 81, 105021 

Kurkela, A., Romatschke, P., Vuorinen, A., & Wu, B. 2010. larXiv: 1006.40621 

Lenaghan, J. T., Rischke, D. H., & Schaffner-Bielich, J. 2000, Phys. Rev., D62, 085008 

Mao, G.-J., Iwamoto, A., & Li, Z.-X. 2003, Chin. J. Astron. Astrophys., 3, 359 

Menezes, D. P., Pinto, M. B., Avancini, S. S., Perez Martinez, A., & Providencia, C. 2009, 
Phys. Rev. C, 79, 035807 



- 19 - 

Menezes, D. P., Pinto, M. B., Avancini, S. S., & Providencia, C. 2009, Phys. Rev. C, 80, 
065805 

Mizher, A. J., Chernodub, M. N., & Fraga, E. S. 2010, Phys. Rev. D, 82, 105016 
Oppenheimer, J., & Volkoff, G. 1939, Phys.Rev., 55, 374 

Orsaria, M., Ranea-Sandoval, I. F., & Vucetich, H. 2011, Astrophys. J., 734, 41 
Paczynski, B. 1992, Acta Astron., 42, 145 

Page, D., Geppert, U., & Weber, F. 2006, Nuclear Physics A, 777, 497 

Page, D., Lattimer, J., Prakash, M., & Steiner, A. W. 2004, The Astrophysical Journal 
Supplement Series, 155, 623 

Papazoglou, P., Schramm, S., Schaffner-Bielich, J., Stoecker, H., & Greiner, W. 1998, Phys. 
Rev., C57, 2576 

Papazoglou, P., et al. 1999, Phys. Rev., C59, 411 

Paulucci, L., Ferrer, E. J., de la Incera, V., & Horvath, J. E. 2011, Phys. Rev., D83, 043009 

Perez Martinez, A., Perez Rojas, H., & Mosquera Cuesta, H. 2008, Int. J. Mod. Phys., D17, 
2107 

Perez Martinez, A., Perez Rojas, H., Mosquera Cuesta, H. J., Boligan, M., & Orsaria, 
M. G. 2005, Int. J. Mod. Phys., D14, 1959 

Perez Martinez, A., Perez Rojas, H., & Mosquera Cuesta, H. J. 2003, European Physical 
Journal C, 29, 111 

Rabhi, A., Pais, H., Panda, P. K., & Providencia, C. 2009, J. Phys., G36, 115204 
Ryu, C.-Y., Cheoun, M.-K., Kajino, T., Maruyama, T., & Mathews, G. J. 2011 



-20 - 

Schwinger, J. S. 1951, Phys. Rev., 82, 664 
Sinha, M., & Mukhopadhyay, B. 2010 
Suh, I.-S., Mathews, G. J., & Weber, F. 2001 
Tolman, R. C. 1939, Phys.Rev., 55, 364 

Weber, F. 1999, Pulsars as astrophysical laboratories for nuclear and particle physics, 1st 
edn. (Bristol: Institute of Physics) 

Weinberg, S. 1990, Phys. Rev. Lett., 65, 1181 

Yakovlev, D. G., Kaminker, A. D., Gnedin, O. Y., & Haensel, P. 2001, Physics Reports, 
354, 1 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



-21 - 




1000 1200 1400 1600 

h, (MeV) 



Fig. 1. — (Color online) Effective magnetic field as a function of baryonic chemical potential 
shown for different central magnetic fields. 
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Fig. 2. — (Color online) Deconfinement order parameter as a function of baryonic chemical 
potential imposing local charge neutrality shown for different central magnetic fields. Dotted 
lines include AMM. 




Fig. 3. — (Color online) Equation of State (pressure of matter as a function of energy density 
of matter) assuming global charge neutrality for different central magnetic fields. Dotted lines 
include AMM. 
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Fig. 4. — (Color online) Particle densities as a function of baryonic chemical potential as- 
suming global charge neutrality for B=0. 
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Fig. 5. — (Color online) Particle densities as a function of baryonic chemical potential as- 
suming global charge neutrality for B c = 5 x 10 5 MeV 2 including AMM. 
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Fig. 6. — (Color online) Adiabatic index as a function of baryonic density for different central 
magnetic fields including AMM. 
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Fig. 7. — (Color online) Moment of inertia as a function of radius for different central 
magnetic fields including AMM. 
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Fig. 8. — (Color online) Mass-radius diagram shown for different central magnetic fields 
including AMM. The dashed and dotted lines have the extra magnetic pressure term of 
B* 2 /247r and B* 2 /8n, respectively, added to the total pressure. 
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Fig. 9. — (Color online) Cooling curves for a 1.1 M© mass star in the presence of differ- 
ent central magnetic fields including AMM. The dashed and dotted lines include the extra 
magnetic pressure term of B* 2 /24tt and B* 2 /8tt, respectively. All curves overlap. 
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Fig. 10. — (Color online) Cooling curves for a 1.4 MO mass star in the presence of differ- 
ent central magnetic fields including AMM. The dashed and dotted lines include the extra 
magnetic pressure term of B* 2 /24tt and B* 2 /8tt, respectively. These two curves overlap. 



31 




Fig. 11. — (Color online) Cooling curves for a 1.93 M© mass star in the presence of different 
central magnetic fields including AMM. The dashed and dotted lines have extra magnetic 
pressure term added to the total pressure of B* 2 /24ti and B* 2 /8tt, respectively. 



